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Within the recently proposed structured FRW model universe the averaged Einstein equations 
are derived. The backreaction turns out to have an interesting behavior. Its equivalent density and 
pressure, being proportional, are negative at early times of the dark ages of the universe, and change 
sign near our present time in our local patch. In addition to explaining the observed dimming of the 
SNIa it leads to new effects for small cosmic redshifts and also to the difference between the local 
and global Hubble parameter. Interpreting the backreaction in the FRW -picture, it is equivalent to 
a time dependent dark energy with w = — 1. 
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Take the universe as it is: homogeneous at large scales 
and inhomogeneous at small scales within structured 
patches. How are the deviations from the standard 
homogeneous cosmic fluid in the matter dominated 
phase of the universe reflected in the cosmological data? 
Recent observational data on SNIa imply a larger dis- 
tance to supernovae than predicted by the conventional 
FRW universe QSHj leading to the term acceleration 
of the universe)!], 0|, and to the concept of the dark 
energy. 

We have recently proposed 6] the Structured 
FRW(SFRW) model of the universe as a first step 
to incorporate the local inhomogeneity of the cosmic 
fluid into a model universe in accordance with the 
observational needs. In the SFRW model of the universe 
the local patches, grown out of the primordial pertur- 
bations, and their backreaction on the homogeneous 
background are modelled exactly as a truncated flat 
Lomaitre-Tolman-Bondi (LTB) manifold embedded in a 
FRW universe from which a sphere of the same extent 
as the LTB patch is removed. As a result of the junction 
conditions the mean density of any such inhomogeneous 
patch, with over- and under-dense regions, has to be 
equal to the density of the FRW bulk[y,0|. Therefore, 
the Copernican principle is in no way violated and we 
are led to a model universe where the local patches are 
distributed homogeneously in the bulk having the same 
mass as a local FRW patch would have, accounting for 
all the structures we see grown out of the primordial 
perturbations within a FRW universe. The analysis of 
the luminosity distance relation in our structured FRW 
model showed explicitly a dimming of objects within 
a patch relative to what it would be inferred from a 
standard FRW universe 0. 

The local inhomogeneous matter dominated patch has 
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a geometrical domain denoted by D and a hypersurface 
£ as its boundary to the FRW bulk. Our calculation 
is based on an exact general relativistic formulation of 
gluing manifolds. The inhomogeneous patch contain- 
ing dust matter is represented by a flat LTB metric 
embedded in a pressure-free FRW background universe 
with the uniform density pt,. We choose the general 
LTB metric to be written in the synchronous comoving 
coordinates in the form0: 
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+ rW], (1) 



where the familiar LTB metric function is now defined 
as R(r,t) — a(t, r).7\The similarity to the Robertson- 
Walker metric as now obvious. The overdot and prime 
denote partial differentiation with respect to t and r, re- 
spectively, and fc(r) is an arbitrary real function such that 
k(r)r 2 < +1 playing the role of the curvature scalar k in 
the FRW universe. Hence, the flat LTB is defined by the 
vanishing of k(r). For a homogeneous universe, a and 
k don't depend on r and we get the familiar Robertson- 
Walker metric. In our SFRW universe, the metric outside 
the inhomogeneous patch, is Robertson- Walker again. 
The corresponding field equations and the solution for 
the flat case k(r) = 0, can be written in the following 
familiar form: 
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where we have introduced g = &M Jf ^ . 
defined as 
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The mass M(r) is 



M(r) = 



2 1 



p(r,t)R 2 dR= ^p{r,t)R 3 , (5) 



where p(r, t) is the density and p, as a function of r and t, 
is an average density up to the radius R(r,t). Note that 
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the volume element in the integral above is not in general 
equal to the proper volume element of the metric, except 
for the flat case k(r) = we will consider 0- In contrast, 
the average density in the patch defined by using the 
proper volume element will be different from the average 
above, except for the flat LTB case. The field equations 
(2-5) are very similar to the familiar Friedmann equa- 
tions, except for the r-dependence of the different quan- 
tities. Furthermore, we assume R'(r,t) = a + ra' > 
to avoid shell crossing of dust matter during their ra- 
dial motion. t n (r) is an arbitrary function of r appear- 
ing as an integration 'constant'. This arbitrary func- 
tion has puzzled different authors who give it the name 
of 'bang time function' corresponding to the big bang 
smg ularityjl GH E3|- 

It has, however, a simple astro- 
physical meaning within our structured FRW universe. 
As R(r, t) is playing the role of the radius of our local 
patch, the time t = t n , leading to R = 0, means the 
time of the onset of the mass condensation or nucleation 
within the homogeneous cosmic fluid. That is why we 
have preferred to use the subscript n for it indicating the 
time of nucleation. As was pointed out in 6], for a realis- 
tic density profile, t n is a decreasing function of the coor- 
dinate r having a maximum at r = 0, i.e. at the center of 
our patch. This means, contrary to the usual interpreta- 
tion in the literature, that t > t n for all < r < ry, = L. 
Therefore, for all times after the onset of mass conden- 
sation within our patch R(t, r) is non-vanishing and for 
times t < t n (r = 0) we have the full FRW without any 
structure. 

Now, without going into the detailed discussion (see 
we know already that, assuming there is no thin 
shell at the boundary of the matching, we must have 

P = Pb, (6) 

where = means the quantities are to be taken at the 
boundary to the FRW bulk. We, therefore, are left with 
the only case imposed by the dynamics of the Einstein 
equations in which the mean density of a local patch is 
exactly equal to the density of the background FRW uni- 
verse: a desired exact dynamical result reflecting the va- 
lidity of the cosmological principle at large, meaning each 
nucleated patch within the FRW universe have the same 
average mass density as the bulk. The total mass in a 
local patch, being equal to the background density times 
the volume of the patch, is distributed individually due to 
its self-gravity, leading to overdense structures and voids 
to compensate it. Assuming again the matter inside each 
patch to be smoothed out in the form of an inhomoge- 
neous cosmic fluid, we expect it to be overdense at the 
center, decreasing smoothly to an underdense compensa- 
tion region, a void, up to the point of matching to the 
background. 

The density distribution within a patch must be such 
that the overdensities of structures are compensated by 
voids. The nucleation time signals the onset of conden- 
sation in the patch which- at least partially- opposes the 



overall expansion. The running of the function t n is cru- 
cial for the expansion history of the patch and therefore 
will influence the luminosity of the structures growing 
within the patch. So far it was shown that t' n < 00. Of 
course, the nucleation time function is related to the ac- 
tual mass distribution for which, taking into account the 
fine structure of the patch including the substructures, 
we have to rely on the overall observations and the mat- 
ter power spectrum [T^ . \l$L ITiL [Tl j . 

We envisage now an averaging process in which the inho- 
mogeneities within the local patch are smoothed out and 
we have again a FRW-type homogeneous modeling of our 
local patch. The traditional way of doing cosmology is to 
take the average of the matter distribution in the universe 
and write down the Einstein equations for it, adding some 
symmetry requirement. One then solves the equations 
G>„ = assuming homogeneity and isotropy of the 

mass distribution as the underlying symmetry. This is 
based on the simplicity principle much used in theoret- 
ical physics. As far as the precision of the observations 
allow, we may go ahead with this simplification. The 
more exact equation, however, is (G M „) = (!),„}. Call- 
ing the difference G M „ — (G M „) = Q^, one may write 
the correct equation as G MJ/ = (T^ v ) + Q^ v . The back- 
reaction term Q has so far been neglected in cosmology 
because of its smallness. Now that measuring Q is within 
the range of observational capabilities we have to take it 
into account. There is, therefore, no need yet to change 
the underlying general relativity or introduce any mys- 
terious dark energy to mimic Q. Of course, the aver- 
aging process is neither trivial nor unambiguous, but it 
is the art of physics to master it. Fortunately, there is 
an averaging formalism, developed mainly by Thomas 
Buchert|16l nA llSl which can easily be adapted to 
our LTB patch, having the same mass as the the FRW 
sphere cut out of it. In this formalism the space- average 
of any function f(t, r) is defined by 




where dV is the proper volume element of the 3- 
dimensional domain D of the patch we are considering 
and Vh is its volume. It has been shown [l^ |20| that in 
such a mass preserving patch the space-volume average 
of any function f(r, t) does not commute with its time 
derivative: 

(fy - U) = (.m -(.nm, (8) 

where the expansion scalar 8, being equal to the minus 
of the trace of the second fundamental form of the hy- 
persurface t — const., is now a function of r and t. The 
right hand side trivially vanishes for a FRW universe be- 
cause of the homogeneity. This fact has far-reaching 
consequences for observational cosmology in our non- 
homogeneous neighborhood. The variation of the Hub- 
ble function with respect to the red-shift is not so simple 
any more as in the simple case of FRW universe 6] . This 
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affects a lot of observational data processing which so 
far has been done assuming homogeneity of the universe. 
Depending on the smoothing width Az, the bins, and the 
matter power spectrum there may be huge effects due to 
the non-commutativity of the averaging process |2 1| . 

The averaged scale factor is defined using the volume 

i 

of o ur p atch D by au = V(t)jj. Now it can be shown 
that El Hi 
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V 



3^ = 3H D . 

a D 
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where we have used the notation &£> = ^ao, and de- 
noted the average Hubble function as Hd- Averaging 
over the local patch means we are taking it as an effective 
FRW patch. Therefore all the derived quantities should 
be based on the average value ap. This is why we take 
the above definition for the mean Hubble parameter and 
not (-), which is different from — . A similar difference 
holds for the second derivative of a: 



a (a) do 



(10) 



Therefore, the definition of the averaged deceleration pa- 
rameter is not without ambiguity, specially because there 
is no nice relation like (9) for the deceleration parame- 
ter. To choose the most appropriate definition, we make 
recourse to the fact that in the averaging process we 
are taking our patch to be homogeneous and FRW-like. 
Therefore, in averaging the redshift as a function a, we 
always encounter ap and its time derivatives ho and ajj. 
This justifies the above definition of the mean Hubble 
parameter and motivates us to make the following defi- 
nition for the deceleration parameter: 
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as was done in the literature so far@, 0, |2^, |2^, |24| . 

Now, we are ready to take the average of the Einstein 
equations in our local patch to see how the mean field 
equations will look like and what are the differences to 
the simple FRW field equations. The emergence of a 
crucial term in the mean Einstein equations, the so called 
backreaction term, is interesting. Buchert's backreaction 
term is defined bv [T^.IT9| 



Q 



((d-(6)f) 



K}-o[r> 



(12) 
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where a is the shear scalar and 8 is the expansion. Al- 
though 9p and Hd are proportional, (8 2 ) and (H 2 ) are 
not. Hence, the relations (13. 14) can not be written 
in terms of H, as was done in|2J|. The averages of the 
Einstein equations using the Hamiltonian constraint and 
the Raychaudhuri equation, taking into account the sub- 
tleties of the observation just mentioned, is then written 



in the following form|16l Il9|: 



(14) 
(15) 



where we have set (p) — pb, the density of the background 
FRW universe, as a result of the junction conditions re- 
flected in the eq.(6). Note that in the so-called Fried- 
mann equation (15) the local Hubble parameter enters 
instead of the global one Hb- The effect of the backre- 
action within the local patch is realized as an effective 
perfect fluid with the equation of state 



PQ = gPQ- 



(16) 



The backreaction term Q can not yet be considered as 
representing dark energy in the FWN-picture. The pb 
appearing in the field equation of our SFRW model is 
the total background energy density, i.e. pb = pu + A, 
where we have chosen pm for the matter density. The 
curvature term is set equal zero to have a simple flat 
universe. But let us first investigate the sign of Q which is 
crucial for the interpretation of these averaged equations. 
As the running of the density and the nucleation time 
t n influence the mean values of the Hubble parameter 
and the shear scalar, the sign of Q is determined by the 
balance between the mean values of the shear and the 
term related to the mean values of the Hubble parameter 
and the expansion scalar in a complex manner depending 
of the running of the density and the nucleation time. 
Given this complex behavior of the backreaction term, 
let us approximate t n in the following way: 



tn - —r 



(17) 



For j3 > the above expansion satisfies all the necessary 
conditions to be fullfilled by t n within the SFRW model 
universe 0. It happens that this behavior corresponds to 
the special parabolic case of ||| formulated in their step 
5 and illustrated in their figure 7. a. The calculation of 
different terms in Q is messy but exactly doable. The 
result is 
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To understand its behavior, we determine its sign for 
two limiting cases: at the onset of nucleation, i.e. t — 
to <C f3L 2 , where the coordinate dependence of different 
quantities has the biggest effect, and at the present time 
t — to ^> (3L 2 . From the exact result of Q we obtain 

t-t ^/3L 2 : Q ~ -50.26^ < 0, (19) 

t - t > I3L 2 : Q ~ +49.75-^ > 0. (20) 

pL z t 

We, therefore, conclude that the backreaction has its 
strongest effect on the onset of mass condensation at the 
beginning of the cosmic dark ages and after the density 
contrast increases to a proper value, where its effective 
density and pressure are negative. It then changes some- 
where, probably at the end of the dark ages, the sign and 
behaves as a normal fluid. 

Negative values of Q at the early stages of mass con- 
densation is a novel effect, even though it is just for a 
local patch. It has the effect of reducing the Hubble pa- 
rameter and producing a negative pressure for a range of 
expansion time with the effect of dimming of the cosmic 
objects in our vicinity. Sometime near our present cosmic 
time the effect reverses and has to lead to new effects for 
small rcdshifts, which are somehow opposed to the accel- 
eration of the universe. Remember that these effects are 
just within our local patch and in interpreting data along 
the light cone one must be cautious. The mere fact that 
local Hubble parameter, Hu, may be less or greater than 
the global one, depending on the behavior of Q, is inter- 
esting and should be taken into account in announcing 
the H values. 

A comparison to the dark energy concept is possible if we 
switch from the SFRW-picture to the FRW-picture. Let 



us denote — Q = A, where A is now a function of time 
and space being most of the time positive. Then we have 
the following equations at our disposal 

Hi = (^) 2 = \ PM (21) 
a D ' 3 

H l = 2 = \(PM+V (22) 

— = -kp M -3A). (23) 

Observational cosmologist, using the FRW-picture are 
used to Hb equation, but taking both the values of 
Hf, and pb = pu + A from observational data. For 
the interpretation of the dimming of cosmic objects, 
however, there is no other way than to use the third 
equation above. In this picture our backreaction can be 
interpreted as a time dependent cosmological constant 
having w = — 1. The actual SFRW picture, however, 
gives us a much wider spectrum of information we 
should be aware of. In the same picture we may say 
that the backreaction produces a dissipative pressure or 
anti-frictional force Q = —A along the line of reasoning 
in|2J|. A realistic structured FRW universe not only 
explains the dimming of cosmological objects but also 
leads to new effects which should be looked for in the 
huge data already existent. 
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